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Abstract 

We investigate the energy evolution of the azimuthal spin asymmetries in semi-inclusive hadron 

production in deep inelastic scattering (SIDIS) and Drell-Yan lepton pair production in pp collisions. 

The scale dependence is evaluated by applying an approximate solution to the Collins-Soper- 

Sterman (CSS) evolution equation at one-loop order which is adequate for moderate Q 2 variations. 

This describes well the unpolarized cross sections for SIDIS and Drell-Yan process in the Q 2 range 

£>y of 2.4-100GeV 2 . A combined analysis of the Sivers asymmetries in SIDIS from HERMES and 

COMPASS experiments, and the predictions for the Drell-Yan process at RHIC at \/~S = 200GeV 

are presented. We further extend to the Collins asymmetries and find, for the first time, a consistent 

description for HERMES/COMPASS and BELLE experiments with the evolution effects. We 

emphasize an important test of the evolution effects by studying di-hadron azimuthal asymmetry 

in e + e~ annihilation at moderate energy range, such as at BEPC at V5 = 4.6GeV. 
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Introduction. Transverse spin azimuthal angular asymmetries in hadronic processes have 
attracted great attentions in recent years. This is not only because the associated observables 
are keen to provide important information on nontrivial hadronic structures, but also because 
they are sensitive to the strong interaction dynamics (l|. The latter involves core feature 
of quantum chromodynamics (QCD): the factorization and universality of the associated 
parton distributions and fragmentation functions, and the energy evolution in hard scattering 
processes. 

Among of these observables the major focuses are the Sivers and Collins asymmetries in 
semi-inclusive hadron production in deep inelastic scattering (SIDIS) and Drell-Yan lepton 
pair production in pp collisions, and di-hadron production in e + e~ annihilation processes. 
The Sivers effects come from the asymmetric transverse momentum dependent (TMD) par- 
ton distribution in nucleon which correlates with the transverse polarization vector S±, 
whereas the Collins effects come from the similar correlation in the fragmentation pro- 
cess associated with the quark polarization. However, these two functions have different 
universality properties: Sivers function differs by sign between SIDIS and Drell-Yan pro- 
cesses [2, [3j, while the Collins function is universal between SIDIS and e + e~ processes. Both 
asymmetries have been observed in SIDIS from HERMES, COMPASS, and JLab Hall A 
experiments [4|-[7|]. In addition, Collins asymmetry has been observed in e + e~ process by 
BELLE collaboration @]. 

The experimental test of the above universality, in particular, for the Sivers asymmetries 
between Drell-Yan and SIDIS, is one of top questions in hadronic physics. However, the 
Sivers asymmetries were observed in SIDIS with Q 2 around 3GeV 2 , whereas the Drell-Yan 
processes will be measured in the range that is greater than 20GeV 2 . In order to consoli- 
date the universality test, the Q 2 dependence of the Sivers asymmetry must be understood 
correctly. The theoretical framework to study the energy evolution of these observables 
has been well developed, where the Collins-Soper-Sterman (CSS) equation I9L 10(] for both 



spin-average and single-spin dependent cross section has been derived |llHl6l|. The CSS 
formalism has been applied successfully to describe low transverse momentum distribution 
of vector boson (Drell-Yan, W/Z) production in unpolarized pp collisions (see, for example, 



Ref. [17]). Early estimate for Q 2 dependence of the SSA [llj was limited to high Q 2 range. 
A recent calculation found a surprising strong evolution effects from HERMES/COMPASS 
energies to typical Drell-Yan energy [18J. This evolution formalism was later applied in a fit 
to HERMES/COMPASS data (3J. The result of Ref. has raised great concerns in the 
experimental proposals, since the predicted asymmetries for Drell-Yan processes would be 
too small due to the evolution. In this paper, we will examine these studies, and carefully 
investigate the Q 2 evolution of both spin average and single-spin dependent cross sections. 



By doing so, we find that the previous study of Ref. [18] over-estimated the evolution effects. 
In particular, the transverse momentum spectrum of the Drell-Yan process in the relevant 
Q 2 range can not be described by the TMD quark distributions proposed in Ref. 15|, [l8j 
(see Fig. 1 below). 

In our calculation, we take an alternative approach, following the original suggestion of 
Ref. 12|, ll3(, by directly applying the CSS equation at one-loop order from low to high ener- 



gies. The one-loop evolution kernel contains a term which predicts a Pt broadening effects 
at higher Q 2 . We will show that this can describe the transverse momentum distribution for 
both SIDIS and Drell-Yan processes, which cover Q 2 in the range of 2.4-100GeV 2 . We extend 
the evolution to the Sivers asymmetries in these processes, and perform a combined fit to the 
HERMES and COMPASS data. The predictions for the SSA in Drell-Yan process at RHIC 



will be updated with the evolution effects. Finally, we will apply the evolution equation to 
the Collins asymmetries in SIDIS and di-hadron production in e + e~ annihilation. 

Collins- Soper-Sterman Evolution. We take the SIDIS as an example, where e(£)+p(P) — > 
e(£')+h(Ph)+X, which proceeds through exchange of a virtual photon with momentum g M = 
£fj, — £'n and invariant mass Q 2 = —q 2 . When Ph± <C Q, the transverse-momentum-dependent 
factorization formalism applies, according which the differential SIDIS cross section can be 
written as 



da(S±) 

= a x 



dxBdydzhd 2 P] 



h±_ 



F uv + e aP SfF. 
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where <To = ^a 2 ra S ep /Q A x (l—y+y 2 /2)xB, and y, Xb, and zu are usual kinematics for SIDIS. 
We only keep the terms we are interested in: Fuu corresponds to the unpolarized cross sec- 
tion, and -Fevers to the Sivers function contribution to the single-transverse-spin asymmetry. 
Fuu and i^ivers depend on the kinematical variables, xb, Za, Q 2 , and Ph_i_, can be written 
into a factorization form with TMD quark distribution and fragmentation functions and 
soft and hard factors. The Q 2 dependence of -FW sivers can be calculated from perturbative 
QCD, and is controlled by the CSS evolution equation, which is easily formulated in the 
impact parameter space. For example, for F? veTS (Q; Phi.) = J 7^2e lPh± ' b ^ Zh F^ vers (Q; b) 1 , we 
have [To!, 

F£ re jQ; b) = F» vcrs (Q ; b)e- s ^ Q > Q ^ . (2) 

The perturbative calculable evolution effect has been included in the Sudakov form factor 
Ssud- In the complete CSS resummation, Q was set at 1/6, and the b* prescription was intro- 
duced: 6* = 6/a/1 + b 2 jb?~ x to deal with the Landau pole singularity |10t j . This necessarily 
introduces a non-perturbative form factor [l0|, which can be determined by comparison to 
the experimental data [17J. 

Alternatively, it was argued in Refs. [l2|, |l3| that we can avoid the Landau pole singularity 
by a direct integration from low to high energy scale, 



Ssud — 2Cf 



Q At, 



dn a s (id) 



In.gUln^ 2 3 



W eg 2 



(3) 



where Co = 2e~ 7B with the complete one-loop coefficients from a recent calculation (20| , and 
both Q and Qo are chosen in the perturbative region. Because of the residual log term in the 
integral, the above Sudakov form factor is not the complete solution to the CSS evolution. 
But, it is a reasonable approximation in the moderate Q and Qq range, in particular, between 
the HERMES/COMPASS and typical Drell-Yan energy regions. For very large Q, such as 
W/Z boson production in pp collisions, we have to take into account higher order corrections 
and back to the CSS formalism. We notice that the same Sudakov form factor applies to 
both Fuu and -F S i vers since they share the same evolution kernel and hard factors in the TMD 
factorization. It works for the Drell-Yan lepton pair production in pp collisions as well 2 . 

Before we study the energy evolution of the SSA, we shall check the above equation 
can describe the spin-average cross sections in the relevant energy range. The majority of 
the SIDIS data from HERMES and COMPASS are in a relative low Q 2 range. Therefore, 



1 Here we only keep the dominant term at low Ph± region, and neglect higher power correction of Ph±/Q- 

2 The difference in hard factors in the TMD factorization does not affect the evolution equation. 
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FIG. 1: Comparison between the theory predictions with experimental data of the low transverse 
momentum distribution of SIDIS at Q 2 = 2.4GeV 2 Ref. 24] and Drell-Yan lepton pair production 
in pp collisions with various Q 2 range Ref. 25]. The scale evolution comes from the Sudakov form 
factor of Eq. ([3]). The predictions calculated from the TMD quark distributions of Ref. [15j] with 
bmax = 0.5GeF -1 and #2 = O.QbGeV 2 are also shown as red curves in the right panel. 



we set the lower scale Qq around these experiments, where a Gaussian assumption for the 
TMD quark distribution and fragmentation functions can well describe the data 2l|, |22 . 
Translating this into the impact parameter space, we parameterize Fuv as 



Fuu(Q ,b) = J2 e\f q {x B )D q 



Sh)e 



-gob 2 -g h b 2 /zl 



(4) 



at Ql = 2AGeV 2 , where / 9 (xb) and D q (z h ) represent the quark distribution and fragmen- 
tation functions following the CTEQ and DSS set Ref 23j at lower scale Q 2 , = 2AGeV 2 . In 
the above equation, go and gh represent the transverse momentum dependence coming from 
the distribution and fragmentation, respectively. In the left panel of Fig. 1, we compare the 
above prediction to the multiplicity distribution in SIDIS from HERMES experiment [24| , 
where we have chosen g = 0.097 and g^ = 0.045. These parameters agree well with those 
in Ref. (H, M ■ 

We can study the Q 2 evolution by comparing to the fixed target Drell-Yan process, 
with Q 2 range from 20 to 100 GeV 2 . To calculate the transverse momentum spectrum for 
this process, we apply the universality of the TMD quark distributions, and the evolution 
equation from Q scale to higher Q. We plot the comparisons between the theory calculations 
with the experimental data in the right panel of Fig. 1. The broadening effects for the Drell- 
Yan processes are well reproduced by the evolution effects of Eqs. (12|3|) . For comparison, we 
also plot the predictions from the TMD quark distributions calculated from Ref. [18j] with 
their evolution effects. Clearly, Ref. [18( over-estimate the broadening effects. It is caused by 
a modification of the non-perturbative form factors used in Ref. [17J in order to describe the 
current SIDIS data, which unfortunately breaks the original predictions for the Drell-Yan 
processes in Ref. (17|. 

Sivers Asymmetries in SIDIS and Drell-Yan. Now, we turn to the Sivers single spin 
asymmetries in SIDIS and Drell-Yan processes. Similar to the above, we parameterize -F s " vers 
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FIG. 2: Theory fit to the experimental data on Sivers single spin asymmetries in SIDIS , as 
functions of x B - left panel from COMPASS [3] and right from HERMES 0|. Q 2 evolution has 
been taken into account from Eq. ([3]). 



at low energy scale Ql = 2AGeV 2 , 
~ b a M 



J2 4 A/™(x) D q {z)e- {go - g ° )b2 - 9hb2/z Z 



(5) 



where M = lGeV is a normalization scale, and we have chosen an additional parameter g s 
for transverse momentum dependence and the fragmentation part remains the same. The 
function Af q (x) = N q x a "(l — x) 13 " 1 q ^ I q f q (x) parameterize the x-dependence of the 

quark Sivers function similar to that in Ref. [2l|. We have the following free parameters: 
g s , a q , (3 q and N q for valence up, down, and sea quarks. Since the data are not sufficient to 
differentiate g s for different flavors, we choose the same g s . We further assume the same sea 
quark parameterization for up, down and strange quarks. 



TABLE I: Parameters {a®} describing our optimum A/j in Eq. (|5|) at the input scale Q 2 = 2.4 GeV. 



flavor i 



N; 



O; 



Pi 



9s (GeV 2 



u 

d 

sea 



0.1±0.004 
-0.5±0.010 
0.5±0.133 



0.88±0.088 
0.23±0.100 
0.43±0.107 



4.1±0.58 

5.3±1.13 

0.016±0.17 



0.024±0.005 
0.024±0.005 
0.024±0.005 



With the above parameterization and the energy evolution effects taken by Eqs. fT2|3|) for 
both spin-average and single-spin-dependent cross sections, we perform a combined fit to 
the Sivers asymmetries from HERMES and COMPASS experiments which scans Q 2 ~2.4- 
lOGeV 2 . We have total of 420 data points, with a minimum \ 2 fit- The best fit results into 
X 2 /d.o.f = 1.67 and the parameters listed in Table I. As an example, we show in Fig. 2 the 
comparisons between the theory calculations and the experimental data as functions of xb 
for COMPASS and HERMES experiments, which demonstrate a consistent description of 
both data. 

Having constrained quark Sivers functions from HERMES/COMPASS experiments, we 
will be able to make predictions for the SSAs in the Drell-Yan processes with the evolution 
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FIG. 3: Predictions for the Sivers single spin asymmetries of Drell-Yan lepton pair production at 
RHIC, v5 = 200GeV, as functions of rapidity for two different mass ranges. As a comparison, we 
also show the prediction without the evolution effects for Q = 4GeV case as dotted line. 
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FIG. 4: The Collins asymmetries in di-hadron azimuthal angular distributions in e + e annihila- 
tion processes: fit to the BELLE experiment at vS = 10.6GeV Ref. [8|], and predictions for the 
experiment at BEPC at \/S = 4.6GeV. 



effects. In Fig. 3, we show that for RHIC experiment at VS = 200GeV, as function of rapid- 
ity with P± integrated up to 2GeV. We have flipped the sign for the quark Sivers function 
because of the nontrivial universality property for the Sivers function. For comparison, we 
have also plotted the prediction without the evolution effects by setting Ss u d = in Eq. (JBJ). 
From this, we see that the evolution reduces the asymmetry by about a factor of 2. This is 
different from that in Ref [18| , where an order of magnitude reduction was indicated for the 
typical Drell-Yan experiments. 

We have done a number of cross checks for the above evolution effects. First, we can tune 
the parameter in the calculations of Ref. 18j to reproduce the P± spectrum of the Drell-Yan 
data, which leads to a much smaller gi = 0.09. With that change, we can describe both 
SIDIS and Drell-Yan data in Fig. 1, and the predicted SSA would be in the similar range 
as ours in Fig. 3. Second, we determine the transverse momentum moment of the quark 
Sivers function (Qiu-Sterman matrix element) from the fit in Fig. 2, and calculate the SSA 



in Drell-Yan process by using the resummation formula in Ref. [20|, neglecting the scale 



dependence of the integrate parton distributions and correlation functions and assuming the 
non-perturbative form factor from Ref. 26| for Drell-Yan process with a mild x- dependence. 
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Again, we obtain the prediction in a similar range as that in Fig. 3. In particular, this 
method provides an important step to matching the SIDIS to Drell-Yan and W/Z boson 
productions in pp collisions. 

Finally, we turn to the energy evolution of the Collins asymmetries. We perform a 
combined analysis to the Collins asymmetries in SIDIS from HERMES and COMPASS 
experiments and the di-hadron azimuthal angular asymmetries in e + e~ annihilation from 
BELLE experiment J8|. Again, we parameterize the Collins function at low energy scale 

Ql = 2AGeV 2 as H^(z,b ± ) = (=%&.) e'^^ 2 /* 2 N q z a «(l - z)^ ^f^ D q {z) and 

\ / a q Pq 

take the evolution effects of Eqs. ([3]). The overall fit is very good, as we show the comparison 
between the theory predictions and the BELLE data. The fitting parameters for the Collins 
functions are listed in Table II, with a x 2 /d.o.f = 1.2. 



TABLE II: 


The fitting parameters 


{«?} 


for the Collins function at the input scale Q 


> 2 = 2.4GeV. 


flavor i 


N- t 




Oii 


A 


9c (GeV 2 ) 


u 

d 


0.34±0.006 
-0.34±0.011 




3.9±0.60 
0.4±0.26 


0.85±0.26 
0.26±0.33 


0.012±0.002 
0.012±0.002 



An important feature we found in this fit is that both favored and disfavored Collins 
functions saturates the positivity bounds. Therefore, it is very important to check the energy 
dependence in other experiments. One idea place is the planed electron- ion colliders |l|, 
where SIDIS processes with wide Q 2 coverage are the major focuses in the proposal. Another 
place is the e + e~ annihilation process at different energies. We suggest to investigate the 
di-hadron azimuthal correlation in e + e~ annihilation at the BEPC of IHEP, Beijing, which 
can reach to the center of mass energy around \/S = 4.6GeV. We show the prediction for 
that energy in right panel of Fig. 3. Earlier experiments at SLAC around the similar ener^ 
range demonstrated applicability of perturbative QCD description of the jet structure |2j 
which shall support to pursue similar studies at BEPC including the Collins asymmetries. 
The initial state radiation events at BELLE can also be used to study this asymmetry in 
various energies [28|. 

Summary. In this paper, we have investigated the energy scale dependence of the spin and 
azimuthal angular asymmetries in hard scattering processes. We applied the Collins-Soper- 
Sterman evolution at one-loop order in the moderate energy range, which can well describe 
the transverse momentum spectrum in existing SIDIS and Drell-Yan data. We focused on the 
energy dependence of the Sivers and Collins asymmetries in these processes, and performed a 
combined analysis with the existing experimental data, including HERMES, COMPASS, and 
BELLE experiments. The non-perturabtive TMD Sivers function and Collins fragmentation 
functions are determined. The predictions for future experiments are also presented. These 
experiments will provide an important test for the TMD universality and strong interaction 
dynamics. We will present more detailed results of this calculation, and the matching from 
SIDIS to Drell-Yan and W/Z boson production in pp collisions in a separate publication. 
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